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jacket structures. One of the most promising of these 
is the tension leg platform (TLP), a semi-submersible 
platform held in place by tension members connected to 
the ocean bottom. 

This thesis evaluates the wave force transfer function 
for the TLP in heave, pitch, and roll motions. This is 
accomplished by assuming the transfer function for a 
Single leg and then extending that result to include the 
effect of phase differences between forces on multiple 


legs, including the effect of wave spreading. 


Cosine squared directional wave spreading is discussed 
and then dismissed in favor of a Gaussian distribution of 
wave direction. Using the Gaussian distribution it is 
shown that when the wave spreading has a standard deviation 
greater than 0.5 radian the exciting force transfer 
functions in heave, pitch, and roll converge rapidly to 
that computed for seas uniformly distributed over all 
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I. “Introduccion 


Present day platforms, in water depths from 10 to 1000 
feet, are primarily of welded steel jacket construction. As 
the present limit of about 1000 feet was approached many 
problems in fabrication and installation were encountered, 
indicating the need for a shift to a new technology. More 
Significantly, with inceasing water depth the natural per- 
iods of jacket platforms increase from about one second in 
100 feet of water to approximately four seconds in 1000 feet 
of water. There is significant wave energy with a period of 
four seconds and dynamic amplification becomes a problem. 
The expected life of the structure is shortened due to high 
cycle - low stress fatigue damage. 

One alternative structure is a tension leg platform 
(TLP). A TLP (see Figure I-1) is a bouyant structure held 
in position by tension members secured to the ocean bottom. 
Because the platform is not rigidly connected to the bottom 
its natural modes of vibration are significantly different 
from a jacket structure of similar height. A TLP in a water 
depth on the order of 1000 feet will have natural periods 
in heave, pitch, and roll of about two seconds, well out of 


the region of high wave energies. 


Since the TLP has only recently been accepted as a 
potential alternative for deep water drilling there are few 
studies of its dynamic response to wave forces. Such preli- 
minary studies are necessary to provide information for the 
design of the operating platforms. Response prediction re- 
quires specification of the modal exciting forces, for example 
in pitch, roll, and heave. For low everyday sea states the 
exciting forces and structure responses can be considered to 
behave in a linear fashion. Linearity of excitation and re- 
sponse suggests that a frequency domain solution technique 
is appropriate; especially when one considers that ocean 
waves are a random process most often described by a wave 
amplitude spectrum. This thesis presents an estimate of the 
modal exciting force spectra in heave, pitch, and roll ina 


way which accounts for structure geometry and wave spreading. 
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Figure I-l Tension Leg Platform 
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II. Problem Statement 


The modal exciting forces on a tension leg platform (TLP) 
will be dependent upon the frequency and the angles of inci- 
dence of the incoming waves. For some angles of incidence 
(at a given frequency) the forces on the various legs will 
act together to produce a larger total force or moment. For 
other angles of incidence the phase differences will cause 
the forces or moments exerted on the various legs to cancel, 
tending to produce a smaller total force. 

In this report the total modal force of the platform in 
heave, pitch, and roll will be divided into two parts. The 
first will be the force exerted upon a single, isolated leg 
by waves of the given frequency. The second will be a geo- 
metric term to take into account the effect of the phase 
differences between forces on each of the legs. This second 
term will be specifically calculated for a square TLP noting 
that the method could be used for any shape. 

The calculations so far imply a unidirectional wave 
amplitude spectrum. Rarely is this observed in practice. 

A reasonable wave spectrum can be approximated as being 


centered about one bearing from the structure and 
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diminishing to zero at bearings 90° either side of the 
central bearing. Two types of wave spreading functions 
will be discussed. The first is the cosine squared spread- 
ing law which, although widely used, allows little if any 
leeway in varying Ene Spreading to fit actual measurements. 
The second is the Gaussian distribution which does allow 
continuous variation of spreading simply by changing the 
standard deviation. 

The modal wave force spectra upon the structure will 
be the integral over all incidence angles of the directional 
wave amplitude spectra, weighted appropriately by the plat- 
form heave force, roll or pitch moment transfer functions. 
The resulting spectra will be normalized in this thesis by 
the results obtained when a uniformly spread wave spectrum 
is used. This is done to demonstrate the effect of wave 
spreading on the heave, pitch, and roll exciting force 


spectra. 


ll 


III. Formulation of Equations 


A. Platform transfer function 

1. Single leg transfer function. The single leg wave 
amplitude to vertical force transfer function, [To(w), is a 
function only of wave frequency. Because the leg has circu- 


lar cross sections there is no dependence on @, the angle of 


incidence. This thesis assumes that [s5(w) is available, 


either experimentally or through another approximation method. 


2. Square platform transfer function. Consider plane 


progressive waves of frequency w coming from an angle g off the 
platform axis. See Figure III-l1. Each leg is thus subject to 
waves of the same frequency but different phase. For deep 


water waves the wavenumber is k = 27 = w*. Thus the wave- 


r g 
length is 


The phase shift in one wavelength is 27. The phase 
shift, die due to a distance, 1, along the axis of prop- 


agation is thus 
2 


2ml; _ w 


Consider the center of the platform as the point of zero 
phase shift. Let d be the leg spacing and a = d//2 be the 
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distance from the center of the platform to the center of any 
leg. 
The distance causing a phase shift for leg 1 is 


l1 =a cos (7 - 8). For leg 2 itis l 


ae ; TO 
1 eo sin(; B). 


T 
Similarly for legs 3 and 4 the lengths are 1, =a sin(7 - B) 
and 1, = -a cos (> - 8). From Equation III-l the phase 
shifts are: 

6 = 


Oa cos(= - 8) 
i J 4 


a) ae 
=-— a sin(= - 8) 
2. 2 4 


2 
w) 
6 = - a sin( 


Tr 
-~ - 8) 
3 g 4 
wy ? T ) 
<a pee a cos (7 - £ i> 7 


Note however that: 


cos (7 - B) cos (7) cos (B) + sin(3Z)sin(8) 


(2/2) {cos(B) + sin(B)} 


sin(? - B) sin(Z) cos (8) = cos (7) sin (8) 


(72/2) {cos(B) - sin(8)} 


Letting T = w*d/g and recalling a = (Y¥2/2)d then 


Equations III-2 become: 
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d = ; {cos(B) + sin(8) } 
" 
T , 
¢ =-= {cos(8) - sin(8)} 
2 2 
jk ; 
¢ = = {cos(8) - sin(§)} 
3 
s =-2 {cos(B) + sin(8)} ae 


3. Heave. The transfer function for heave for the 
total structure, Pee is the sum of the contributions of g 
each leg. Taking phase into consideration each leg contri- 


butes T, exp{i(wt+$,)}. Thus the quantity of interest, the 


square of the magnitude of I,,, is 


lp. /2 2 \r |? Jetoty edn On aides 2. 


ht 


Incorporating the identity er? = cos(4) + i sin(9), 
summing real and imaginary parts, and taking the magnitude 


squared of the exponentials this becomes 


ile = bee ({cos($¢,)+cos($,)+cos(¢;)+cos(¢,) }? 


+{sin($1)+sin($,)+sin(¢,)+sin(¢,) }?) III-4 
Note however from Equation III-3 that 6, = -$, and 
6, = ->,. Thus we have 
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cos(¢,) = cos(¢$,) sin($,) -sin($,) 


cos($,) = cos($,) sin($,) = -sin($,) 


and Equation III-4 becomes 


Was Ir,|* {2 cos(o,) + 2 cos(,) }? 


ITnel? = 4 pals {cos($,) + cos(o,)}?. III-5 


By replacing $, and $, with their values from Equation 


III-3 and setting R = [Moral <// jem this becomes 


4 (cos {= (cos 8 +sin 8 )} + cos{$(cos 8 -sin B )})? 


*h 


w 
i! 


T 
h 4 {cos (> cos B) cos (5 sin 8) - sin(5 cos 8) sin(3 sin 8) 
a cos (5 cos 8) cos (5 sin 8) + sin(5 cos 8) sin(> sin 8) }? 
Ry = 4 {2 sae cos Aces sin 8) }? III-6 
2 2 


A form often more suitable for computer work is obtained 


by applying the identity cos* A = 4 + k cos 2A yielding 


ITnel? = 4 (1+ cos(T cos 8)) (1+ cos(T sin 8)) ie 
III-7 


which can be implemented by recalling that T= w’d/g. 
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4. Horizontal force assumptions. The roll and pitch 


moments experienced by a TLP will be dependent upon horizon- 
tal as well as vertical forces. In this report it has been 
assumed that the total moment producing effect of the hori- 
zontal force is small compared to the vertical forces. This 
is perhaps not always the case but is assumed in this discus- 
sion because, while the horizontal wave forces may be large, 
they are distributed along the legs, both above and below the 
center of rotation of the platform. The force on the leg 
near the surface is large but with a small moment arm. The 
forces below the center of rotation will be smaller but will 
extend further down the leg and thus have a large moment arn, 
producing a moment similar in magnitude but opposite in sign, 
and will tend to cancel the other. Thus the moment producing 
effects of the horizontal forces are assumed to be signifi- 
cantly smaller than those of the vertical forces. When the 
location of the center of rotation and the roll and pitch 
natural periods are known, as they would be for any particu- 
lar TLP design, the horizontal contribution to the roll and 
pitch moment could be calculated. It can be shown that the 
phase angles encountered when adding up the moments generated 
by the horizontal forces on each of the four legs would be 


the same as calculated previously for heave. Therefore the 
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results computed for heave will serve as a check on the qual- 
itative behavior of the role of horizontal forces in the 
generation of roll and pitch exciting forces. 

5. Pitch and roll. Pitch and roll for a square plat- 
form are exactly symmetric. Thus only roll will be discussed 
here with the understanding that pitch could be obtained by 
shifting the coordinates by 90 degrees. Taking the direction 
of positive roll moment to be as shown in Figure III-1l the 


total transfer function for roll, squared, would be 


fh 
tl 


[P_|? | (4/2) 08644 (4/2) ei 674 (-a/2) e994 (-4/2) ef 9* |? 


2 ; 
ale ; JetOt4zel2-eids_cids | 2 


Let R_ = IT eel°Z|T|*%. Incorporating the identity 
eld = cos ? + i sin ? , summing the real and imaginary 


parts, and taking the magnitude squared this becomes 


2 


R = 


- {(cos $,+ cos ¢,- cos $,- cos o,)? 


ih 1 Oz 


+(sin $,+ sin $,- sin $,- sin $,)?} III-8 


However, from equations III-3 $6, = -¢, and $3 = -$, . 
Thus cos $, = cos ¢, , cos $,; = cos d2 , Sin >, = -sin d, . 
and sin $3;= -sin $¢, . Therefore the cosine term in 
Equation III-8 reduces to zero and the sine term becomes 


(2 sin $, + 2 sin 9,)7, leaving 
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= 2 ° ° Z 
R, = d° ( sin’ $j+esin- oe 


Incorporating the values from Equations III-3 this becomes 


ae (sin(3 (cos 8B +sin 8 ) }+sin{-3 (cos 8 -sin Bg )})? 


y°) 
uN 


R = @? (sin (2cos 8) cos (ssin B) + cos (cos 8)sin(<sin B) 


~ sin(-5cos 8) cos (sin B) + cos (~Scos 8)sin(=sin B) }* 


4) 
ut 


T 
a? (sin(5cos 8) cos (ssin 8B) + cos (cos 8) sin(fsin B) 
- sin(2cos 8) cos (2sin 8) + cos (zcos 8)sin(Ssin B)}? 
2 270 5 ol. 
R = 4 a* cos (Scos &8)sin (Ssin B) III-9 


By applying the identities cos*A = 5 + 5 cos 2A and 
sin*A = 5 - 5 cos 2A this yields a form more suitable for 


computer work. 


Ir_, |* = a* {1 + cos(T cos 8)}{1 - cos(T sin 6)} be 


EG 
III-10 
As mentioned previously pitch is simply rotated by 90 


degrees. Thus cos 8 becomes sin 8 and sin 8 becomes -cos 8 . 


The transfer function for pitch for a square platform is thus 


lr 


Sal = d* {1 - cos(T cos 8)}{1 + cos(T sin 8)} Waiee 


IIIT-1l 


By) 


iain Wave spectrum. 


The wave spectrum incident on the platform varies with 
both frequency and direction and is denoted S,(,B,B,). 
8, is the center of an arc of 180° from which the waves 
are coming. 

For simplicity of calculations it is often assumed, as 
it shall be in this report, that the total wave spectrum 
can undergo separation of variables. Thus we can express 


the spectrum as 


S_,(@-8.B,) = S , (a) F(6,8,) 


where the entire frequency dependence is contained in 
S,(#) and the entire angular dependence is in F(8,8,). 
This section deals with approximations for the angular 
dependence, F(§8,8,). 

1. Cosine squared wave spectrum. An often used method 
for describing the angular spreading of waves is the cosine 
squared wave spectrum (see Figure III-2). It has the 
advantage that it vanishes for values of 8 equal to plus and 


minus 7/2. The formula for the cosine squared spectrum is 
S,,§@,8,8,) = C p(8,By) Si(w) s -m/2 < B-B, < 1/2 


where 8, is the center of the wave spectrum and C is a constant 


20 


used to normalize to unity the area between -1/2 and 1/2. 


It is derived as follows: 


t/2 
l1 = f C cos*(8) as 

-7/2 
1 t/2 = 
- = 2 f cos*(8) dB = 2 (=) 
Cc 0 4 
c=2, 

T 


Thus the equation for the cosine squared spectrum is 


N 


S_(@8,B,) = = Si(w) cos*(B-By) + -5 < (B-8) < 5 


Lit—12 


4 


This spectrum suffers from one major disadvantage. 
It cannot allow for variations in the directionality of the 
waves. This problem can be eased somewhat by considering 
cosine to the fourth, to the sixth, and so forth but this 
results in only step changes in the directionality. Methods 
have been designed to provide continuous changes but they are 
somewhat cumbersome and comparison is not apparent. 

2. Gaussian wave spectrum. The method used in this re- 
port to achieve continuous changes in directionality is to 
consider the waves to be spread about the central angle 6, Pig 


a Gaussian fashion with standard deviation o. The reader not 


a 


-1/2 0 (B-8,) 1/2 


Figure III-2 Cosine Squared Spreading 


o=0.5 


S ila PN 


-1/2 0 (B-Bo) 1/2 


Figure III-3 Gaussian Spreading 


familiar with Gaussian (or normal) distributions may wish to 
refer to any good probability text for an explanation more 
detailed than the following. 

The probability distribution of a Gaussian process with 


mean 8, and standard deviation o is described by 


p(8,B,) = 





exp{-(8-8,) 7/(207) } 
/27 oO 


A spectrum of waves defined by a Gaussian process would 


thus be 


S,'%,B,Bo) = C p(B,B,) So) % =/2<8-8, <1/2) et Toe 


where again C is a constant to normalize to unity the area 
between -71T/2 and 71/2. It is derived as follows, referencing 


B, to zero: 


1/2 
1= {| c p(8,Be) 4B 
-1/2 


1/2 


== f p(8,Be) a8 
-1/2 


This integral is evaluated using the properties of the 


Gaussian distribution and yields 


ak <a 
2 (55) -1 


x 
= db _¢2 
where 4$(x) = is exp(-t*/2) dt 


Z3 


is known as the normal distribution function and is tabulated 

in probability texts and math tables. Thus C is a function 

of the standard deviation and shall hereafter be denoted C(o). 
Unlike the cosine squared spreading the Gaussian distri- 

bution does not vanish for 8 equal to plus and minus 1/2. 

This is why C(o) varies with o. It must account for the 

cropping of the "tails" of the distribution for |8| > 7/2. 
Some simplifying approximations can be made however. 

As can be seen from Figure III-3 when o reaches a value of 

2 or larger the angular spectrum approaches a constant 

value. That value is 1/m (in order to achieve unity area). 


Thus Equation III-13 becomes 


= 


TT TT 
S,(@,8,B9) a S _,) ; "5 < B=Bige < > EOL {O0>2. 


For o < 1 the area in the tails diminishes rapidly 
leaving almost all of the area between <-1/2 and 7/2. As 
can be seen in Figure III-4 C(o) rapidly approaches unity. 
For o=1.0 C=1.13. For o=0.9 C=1.09. So foro <l 


Equation III-13 can be approximated by 





z 1 =(8-8,)* 
S,(w,B,Bo) rz S fo) = = exp{—s—z-—} 


For values of o between l and 2 the entire exact formula 


must be used 
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C(o) 


Oc25 1.0 ees 220 
o (radians) 


Figure III-4 Correction Factor for Loss of Gaussian Tails 
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2 
1 - } ap( aS) I1lI-14 





s_(w/88,) = S00) GO Vo ore 20 
20 


3. Comparison of cosine squared with Gaussian. Since 


many people are familiar with cosine squared spreading it 

is useful to discuss for what values of o the Gaussian would 
approximate the cosine squared. Although the two functions 
are quite distinct there are several possible methods to 
compare them. 

The normalized cosine squared distribution has a value 
of 2/m = 0.637 for B = B,-. For Gaussian spreading this would 
be C(a) / (Y27 o) at B = Bo. Ifo is assumed less than l 
then C(o) = 1. Equating the two values yeilds o = 0.627. 

When 8B = 8, + 1/4 the cosine squared distribution is 
half of its maximum value. This situation can be matched 
with a Gaussian yielding still another "equivalent" o for 


the cosine squared. 


C(w) prin /ay", ea C(w) 
Won og 8 207 7 : /2n o 





exp{—Z7) = (0.5) 


- (7/4) ? ene 


ott in(0.5) 


0.667 


Q 
i 
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The two values obtained are approximately equal. Thus the 
Gaussian distribution with o=0.65 is very similar to the more 
familiar cosine squared distribution. A smaller o means a more 


unidirectional wave front and a larger o means more spreading. 


C. Platform exciting forces. 


Each differential element of the spectrum S 648 ,B,) ap 
induces mean square force or moment upon the platform through 
the appropriate platform modal force transfer function, 
magnitude squared, IT, (w,8) | *. Thus the total force or 
moment spectrum is given by 

TT 
H(w,Bo) = f |T (w,8) |* S_(w,8,8,) dB III-15 
-T 

It is important to note that this modal force spectrum 
is dependent upon 8,, the center of the directional spectrum. 

1. Square platform, heave force. The heave force 
spectrum for a square platform in normally distributed 
waves is derived by substituting Equations III-7 and III-14 


into Equation III-15 yielding: 


By, +1/2 
H,(@,Bo) = S le 4{l+cos(T cos 8)}{l+tcos(T sin 8) }x 
By-1/2 
AL 1 -(B-Bo) *) ap 
S fo) ‘TeE-1 Vere =} expt — — 
0 


27 


C(o) 
Hy (a8) = 4 ITI" S30) Tay 


x 





ig = (GEA 
if {l+cos(T cos 8)}{l+cos(T sin 8) } oe ee ae ag III-16 
By-1/2 


Setting Tyhlw,8 ) equal to the integral in Equation 
III-16 and expanding the integrand yields 
Beat 2 


I (a,Bo) = f {1 + cos(T cos 8) + cos(T sin B) + 
Bi -m/2 


es pea) 2 
cos(T cos 8)cos(T sin 8) } exp (Eau) dp III-17 


& 
where H,(w,B,) = 4 [pal i S, (w) a I (w,B,) - 





A closed form solution of Equation III-17 does not 
exist. Thus a numerical integration was necessary 
( see Section IV). 

2. Square platform, roll moment. Similarly the roll 
moment spectrum for a square platform in normally distributed 
waves is derived by substituting Equations III-10 and III-14 


into Equation III-1L5 yielding: 
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Batis 2 


H (a ,Bo) = f 
By-1/2 


ad’ {l+cos(T cos 8)}{l-cos(T sin 8) } | Coie 


= (8 =o 
= a dB 


{ I Fe ent 


USE «oem 
20 





H (Bq) d \T Sa) 5a 


es (B-By)? 
f  {l+cos(T cos 8)}{l-cos(T sin 8)} exp(— = dg III-18 


8,-1/2 


Setting I_(w,8,) equal to the integral in Equation III-18 


and expanding the integrand yields 


B,+1/2 
I_(w,85) = f {l+cos(T cos 8)-cos(T sin 8)- 
8,-1/2 


cos(T cos 8)cos(T sin 8)} exp{—E=52 3 ag III-19 


C(o) 
H_(w,8,) = ol ee Eley a I(w,Bo) . 





where 


Again a closed form solution does not exist and the same 


numerical integration technique was applied ( see Section IV) 
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3. Square platform, pitch moment. As mentioned 


previously, for a square TLP the pitch moment is simply 


rotated 90 degrees from the roll moment. Thus 





2 2 C(c) 
H,(o,B,) = d° |T.]* S$, (a) a I ,(a 1B) 


where TI (w,8,) = I (w,B,+m/2) . 


4. Normalization. In order to compare the effect of 
the wave spreading for different configurations a method 
of normalization was developed. The normalization factor 
was taken to be the force or moment frequency spectrum 
resulting from waves of equal amplitude from all directions 
in one half of the plane. Thus S,(w,8,80) equals (1/m) S, (w) 
for (By-1t/2) < B < (8)+n7/2) and zero elsewhere. The 
normalization frequency spectrum becomes 
Botnr/2 
Hy (ar Bo) = f IT (a,B) |? (l/T) S, (o) ag. 
By-T/2 
Because of the symmetry of a square platform and 
the uniformity of this spreading it makes no difference 
where the center of the spectrum is. For convenience 


8, was taken to be Zero and thus the normalization spectrum 


becomes 
l T/2 
Hy fo) = = silo) fs IT, (7,8) |? ap TII-20 
-1™/2 
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which is simply the mean square value of ry times S §) « 


For heave and roll they are as follows 


1/2 

Bayes Ir_|* fs {l+cos(T cos 8)}{l+cos(T sin 8)} dg 
h T ON Ss ~n/2 

III-21 
q? 1/2 

Hp, (w) == S, (w) lr_|* fs {l+cos(T cos 8)}{l-cos(T sin 8)} dp 

Ds Tr ON S 

-7/2 

ERe—=22 

The normalized spectrum for heave is thus 
H (w,Bo) 
ee >| _° = 
HN (w,8,) Ca) litezs 


Expanding the numerator and denominator and cancelling 


the common terms S_(w) and weal this leaves 





an opt (w,B,) x 
HN (a,Bq) oe) III-24 
where 
H, (o,8,) 
T1, (w,B,) = S (a) Ir [= 
n S 
Botn/2 
“= os if {l+cos(T cos 8)} x 
y B,-1/2 
{l+cos(T sin 8)} exp{—'B-Bal } dp 
2) ape as = 
h S, §@) is 
4 n/2 
= (~) f[ {l+cos(T cos 8)}{l+cos(T sin 8)} dap 
-7/2 
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Similarly for roll 


H.(w,8, ) Mm) (0,6. ) 
HN. (w,8o) i Hoy (#) ~ 3 (w) 
r 


where 


Hy (w,Bo) 


I1 _(w,Bq) = s,(@) Il? 


Bo+t/2 
C(a) s {l+cos(T cos 8)} x 
of B,~/2 





Hl 
4 te 


{l-cos(T sin 8) } exp {Epo ap 


H,_. (w) 


13 _(o) = Sta) [Fel® 


a2 1/2 
= f[ {l+cos(T cos 8)}{l-cos(T sin 8)} dg. 
-17/2 
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IV. Numerical Prediction of Exciting Forces 


A. Discussion. 

Because of the complicated forms of the integrands of 
Equations III-17 and III-19 no closed form solutions to 
the integrals were found to exist. Numerical methods were 
employed to determine the variation of resulting forces 
caused by the variation in input parameters. Appendices 
I and II provide program listings and an explanation for 


their use. 


B. Assumed values. 

In order to use numerical methods of calculation 
values must be assumed for the relevant parameters. In 
this thesis the initial values assumed come from two 
Amoco papers published in the “Journal of Petroleum 
Technology." (See Reference 1.) These papers discuss 
a proposed vertically moored platform (VMP) which is 
identical to the concept of a TLP. Amoco's VMP is proposed 
to have leg spacings of 160 feet and is estimated to have 
heave, pitch, and roll natural periods of approximately two 
seconds. Because of the popularity of the cosine squared 


wave spectrum a comparable standard deviation of 0.65 
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radians was chosen. 

The parameter of interest in all of the following 
analyses is HN, the ratio of the magnitude of the transfer 
function for Gaussian wave spreading to the magnitude of 
the transfer function for uniform wave spreading. 

When the Amoco TLP was analyzed in normally distributed 
waves with a 0.65 standard deviation the results were some- 
what surprising. The variation in HN from unity in heave 
was less than one percent for all values of wave spectrum 
incidence, 8,, (see Table IV-1). A much greater variation 
of HN for roll was expected but HN still remained within 
ten percent of unity (see Table IV-2). Clearly for this 
case an initial assumption of uniform wave spreading and 
the resulting simplifications in calculation of the force 
spectrum will provide very useful first pass information 
to the designer. 

However the Amoco proposal is by no means the only 
viable configuration for a TLP. One must consider the 


variation in HN as other parameters vary. 


C. Variation in leg spacing. 


The most obvious variation ina platform is the spacing 


between the legs. With the 160 feet of the Amoco platform 
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as a reference, leg spacings ranging from 100 to 250 feet 

were considered. Analysis results for heave are presented 

in Table IV-3 and for roll (and hence pitch) in Table IVv-4. 
For all of these configurations HN for heave is still 

within one percent of unity. Roll suffers somewhat but 

the calculated HN stays well within 20% of one. For initial 

design work this accuracy is sufficient and the assumption 


of uniform spreading would be justifiable. 


D. Variation in period. 


Because the wave spectrum is continuous in frequency 
waves of all periods will be incident upon the structure. 
However, the period of paramount concern is the natural 
period of the structure. This varies from structure to 
structure and is somewhat controllable. Again taking the 
Amoco value, two seconds, as a reference, periods from one 
to three seconds were considered. The leg spacing was 
assumed to be 160 feet and the standard deviation 0.65 
radians. Analysis results for heave are tabulated in 
Table IV-5. The results for roll however vary enough to 
be visible on a polar plot of HN vs. 8,, the center of 
the wave spectrum, and are shown in Figure IV-6. When 


possible a plot has been preferred to a table because of 
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the greater amount of information compactly presented. 

For heave the variation in HN remains at about one per- 
cent reaching a high of two percent for a three second per- 
iod. For roll HN shows marked fluctuation but again stays 
within twenty percent of unity. While the natural period 
has significant implications with respect to the periods of 
Significant wave energy it does not appear to significantly 
alter the modal force transfer functions as long as some 


directional spreading is present. 


E. Variation in standard diviation 

The final parameter to vary is the standard deviation. 
The standard deviation for uniform wave spreading is infinite 
while that for unidirectional waves is zero. However, as 
shown in Section III-B standard deviations above two radians 
approximate the condition of uniform wave spreading. There- 
fore almost all of the effect of varying standard deviation 
occurs between zero and two. One would thus expect a signi- 
ficant change in HN as standard deviation changes, particu- 
larly for small values. 

This expectation is correct. Standard deviation was 
varied from 0.125 to 2 radians for both heave and roll 
(D = 160 feet, Per = 2.0 seconds). The results are plotted 
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in Figures IV-7 and IV~-8 for heave and roll respectively. 
Consider heave (Figure IV-7). For standard deviations 
of 1.0 or 2.0 radians HN stays within one percent of unity. 
When o = 0.5 the variation becomes five percent. For 
o = 0.25 it exceeds twenty percent and when the standard 
deviation is 0.125 radians even the pattern of the curve 
changes. For roll (Figure IV-8) the results are similar, 
although somewhat more pronounced. 
However, if one can allow a 15% error in an initial de- 
Sign approximation then the uniform wave spreading assumption 
1s valid for wave directional spectra with standard deviations 
greater than approximately 0.5 radian. Note that cosine 


squared spreading is included in this category. 


F. Summar 

Numerous other cases were considered, simultaneously 
varying leg spacing, period, and standard deviation. In all 
cases where o > 0.5 the resulting normalized amplitude- 
platform force transfer function had a magnitude within 
20% of unity. Thus for waves with any reasonable amount of 
spreading sufficient accuracy for initial design calculations 
can be obtained by using the mean square value of [(w,8). 
This result allows for a simpler initial calculation of 
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olatform response to random waves in the manner of 
Vandiver.? For situations when o < 0.5 or when more 
accuracy is required Vandiver's method can still be 
employed using the computer programs of this thesis to 


determine HN (Vandiver's Cl). 
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SQUARE TENSION LEG PLATFORM, HEAVE EXCITATION 


WITH GAUSSIAN WAVE DISTRIBUTION 


LEG SPACING = 160 feet 
PERIOD = 2.0 seconds 
STANDARD DEVIATION = 0.65 radians 


Correction factor due to loss of tails = 1.016 


Factor due to uniform spreading (I3) = 3.94 
Il HN, =I1/13 
3.90 0.990 
391 0.992 
3.92 0.994 
a.95 0.997 
3.94 1.000 
LPs Js! 1.002 
3396 1.004 
a0 1.006 
3.96 1.006 
3.97 1.007 


Table IV-l Normalized Platform Heave Exciting Force 


(D=160 feet, Per=2.0 sec, o=0.65 rad.) 
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B 


10 
15 
20 
Z> 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 


SQUARE TENSION LEG PLATFORM, ROLL EXCITATION 


WITH GAUSSIAN 


LEG SPACING 
PERIOD = 


STANDARD DEVIATION 


Correction factor due to 
Factor due to uniform spreading (I3) = 0.235 E+05 


je 


0.216 
0.216 
OF 74 
0.218 
Oezis 
Onee22 
0.224 
0.227 
Orzo 
0.234 
0.238 
0.241 
0.244 
0.248 
OZ oe 
0.254 
O25 5 
0.256 
0.257 


WAVE DISTRIBUTION 


160 
2.0 
=0.65 


E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 
E+05 


feet 
seconds 
radians 


loss of tails = 1.016 


HN =I1/13 
rE 


0.919 
0.920 
0.922 
05925 
0.933 
0.943 
0.954 
0.967 
0.981 
0.996 
etek 
1.026 
1.040 
i053 
1.068 
1.079 
1.086 
sah Te 
1.092 


Table IV-2 Normalized Platform Roll Exciting Moment 


(D=160 feet, Per=2.0 sec, o=0.65 rad.) 
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SQUARE TENSION LEG PLATFORM, HEAVE EXCITATION 


WITH GAUSSIAN WAVE DISTRIBUTION 


PERIOD = 2.0 seconds 
STANDARD DEVIATION = 0.65 radians 


CORRECTION FACTOR DUE TO LOSS OF TAILS = 1.016 


LEG SPACING HN), HNy, 
(feet) (min. ) (max. ) 
100 0.997 1.003 

150 0.998 1.003 

200 0.994 1.006 

250 0.995 1.004 


Table IV-3 Normalized Platform Heave Exciting Force 


(Various leg spacings, Per=2.0 sec, o=0.65 rad.) 
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SQUARE TENSION LEG PLATFORM, ROLL EXCITATION 


WITH GAUSSIAN WAVE DISTRIBUTION 


PERIOD = 2.0 seconds 
STANDARD DEVIATION = 0.65 radians 


CORRECTION FACTOR DUE TO LOSS OF TAILS = 1.016 


LEG SPACING HNy HNy 
(feet) (min. ) (max. ) 
100 | 0.978 1.027 

150 0.933 1.067 

200 0.882 1.114 

250 0.846 L260 


Table IV-4 Normalized Platform Roll Exciting Moment 


(Various leg spacings, Per=2.0 sec, o=0.65 rad.) 
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SQUARE TENSION LEG PLATFORM, HEAVE EXCITATION 


WITH GAUSSIAN WAVE DISTRIBUTION 


LEG SPACING = 160 feet 
STANDARD DEVIATION = 0.65 radians 


CORRECTION FACTOR DUE TO LOSS OF TAILS = 1.016 


PERIOD HN, HN 

(seconds) (min. ) (max. ) 
1.0 0.997 1.001 
25 0.996 1.005 
2.0 0.990 1.007 
225 0.994 1.004 
3.0 0.980 1.024 


Table IV-5 Normalized Platform Heave Exciting Force 


(D=160, Various periods, o=0.65 rad.) 
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Figure IV-6 Normalized Platform Roll Exciting Moment 


(D=160, various periods, o=0.65 rad.) 


AG 





= Hy (a, By) 
8 hy Se, 


HN 


Figure IV-7 Normalized Platform Heave Exciting Force 


(D=160, Per = 2.0 sec, various o) 
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Figure IV-8 Normalized Platform Roll Exciting Moment 


(D=160, Per = 2.0 sec, various o) 
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Vie Suggestions for Experimental Verification 


A. Single leg. 

Because the results of the computer program provide a 
transfer function between forces on a single leg to platform 
forces and moments, the single leg forces must be obtained. 
This can be done theoretically or experimentally. This 
report will deal with experimental determination. 

The single leg is assumed to be of circular cross sec- 
tion. Thus there is no angular dependence and testing can 
be performed in a ship model towing tank with a wave gener- 
ator. In the proposed test the single leg model will be held 
rigid with load cells measuring the forces as a result of 
incident waves. These forces will be normalized to forces 
due to scaled unit amplitude waves using a linear assumption 
for the force to wave amplitude relation. By performing 
this test with a variety of wave frequencies a plot jongiaG@® 
can be obtained. From this the computer will predict 


platform heave force and roll moment transfer functions. 


B. Square TLP. 


To completely test the results of this study would re- 


quire the ability to generate a directional wave spectrum. 
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This ability is not available. However, the total transfer 
function, Pi fw,8), can be experimentally verified for heave, 
meted, and roll. 

Once the results for a single leg are obtained a 
Similar experiment with the total platform can be performed. 
The model should be held rigid with load cells measuring 
forces in all six degrees of freedom. The direction of 
wave incidence can be varied by turning the model. 

This experiment will generate a family of curves of 
Pi (w,3), one curve for heave and one for roll for each angle 
of incidence 8. These curves can then be compared to those 
predicted using the measured forces for [,(w) and the heave 
and roll algorithms developed in this paper, Equations III-7 


and III-10. 
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VI. Conclusions 


ip Assumptions and approximations. 


1. Linearity. The forces computed in this” papemmane 
those exerted by waves on a motionless structure. The 
forces are assumed to be linearly related to wave amplitude. 
Furthermore the forces exerted on a motionless structure 
can only be useful in predicting dynamic response if the 
equations of motion are linearized. For a TLP the equations 
of motion will be linear for low to moderate sea states. 
This is due to the large dimensions of the major members and 
to the relatively small response amplitudes. Damping terms 
for the TLP may also be linearized because the damping 
ratios are expected to be small in heave, pitch, and roll. 

2. Effects of cross members. The model of a TLP used 
to derive the transfer functions had four legs acted upon by 
wave forces, and nothing else. Real TLP's must have cross 
members between the floatation legs to provide structural 
support. These would also be subject to wave forces. 7 

The diameter of a cross member will be much smaller than 


that of the leg. In the Amoco platform referred to in Refer- 


ence 1 the cross members have an outer diameter of 5.5 feet 
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while the floatation leg's diameter is 30 feet near the sur- 
face where the large wave forces are. This is a factor of 
about 5% difference. The displaced volume is proportional 
to the cross sectional areas which are almost a factor of 30 
different. In addition the support members are mostly near 
the bottom of the legs, away from the largest wave forces. 
The parts that are shallow are close to the legs and would 
react in phase with them. Thus the effect of the waves on 
the cross members can be assumed to be negligible with 
respect to the effect on the legs. 

3. Wave disturbances of the legs. Another aspect not 
considered in this model is the effect of one leg upon the 
waves incident upon the others. In the Amoco example the 
outer edges of the legs at the water surface are 130 feet 
apart which is only about four times their 30 foot diameters. 
There could be some shielding or reflecting effects present, 
Causing the results of this study to be in error. As the 
ratio of leg spacing to leg diameter increases these errors 
will decrease. 

The primary purpose of this thesis was to evaluate the 
influence of wave spreading on the heave, pitch, and roll 
force spectra. The most important geometric effect is the 


relative phases of forces exerted on the legs. In cases 


50 


that the wave spreading is shown to be sufficient to average 
out the relative phase effects, then it can certainly be 
argued that interaction effects between legs would also 


be suppressed. 


B. Utility of results. 
1. Initial design approximation. For wave directional 


Spreading with a standard deviation of greater than 0.5 
radians the heave force or roll or pitch moment imposed 
upon a TLP can be approximated within 20% to be that imposed 
by uniformly spread seas. By incorporating this into the 
method used by Vandiver’ an initial estimate of the response 
of a TLP can be obtained which would be applicable in 
preliminary design. 

2. More accurate approximation. If the standard 
deviation of the wave spectrum is less than 0.5 radian or 
if an answer with more than 20% accuracy is required then 
the computer program outlined in Appendices I and II can 
be used to determine HN (Vandiver's Cl). This can then 
be incorporated into the calculations of the modal response 
of the structure. 

3. Variation to fit other structures. Results for 


other square TLP's can easily be obtained simply by changing 


yal 


the values assigned for the leg spacing and the period of 
interest. If the structure of interest 1s not square the 
same concept can be applied by changing the program defini- 
tions of the transfer functions. In this manner results 
can be obtained for any shape of TLP (square, triangular, 
rectangular, etc.) with any leg spacing and any natural 
period. 

4. Wave directional spectrum variation. The greatest 
advantage of modeling the wave directional spectrum as a 
Gaussian process is the continuous variation in spreading 
it provides. The standard deviation of the wave spectrum 
at a proposed location for a TLP could be measured under 
various conditions, noting that the periods of the signi- 
ficant wave energy would also change. This information, 
standard deviation and period, could then be input to this 
program to estimate platform forces and moments under those 
conditions. The final information of forces, moments, and 
responses would then be applied to sizing the tension 
members for adequate tensile and fatigue strength and to 
ensure the platform motions do not cause equipment malfunc- 


tions or personnel discomfort. 
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Ox 


HN 


163 
oi 


NOMENCLATURE 
distance from center of TLP to center of any leg 


constant to normalize the area of the wave 
directional spectrum 


leg spacing of TLP 
acceleration of gravity 
aS a subscript, denotes heave 


frequency spectrum of the total force or moment 
on the TLP due to wave spreading 


frequency spectrum of the total force or moment 
on the TLP due to uniformly distributed waves 


normalized frequency spectrum of the total force 
or moment on the TLP 


integral portion of Hy dependent on wave spreading 


numerator of HN, dependent on platform geometry 
and wave spreading 


denominator of HN, dependent only on platform 
geometry 


wavenumber 


Gaistance from center of TLP to projection of leg i 
on the direction of wave propagation 


as a subscript, denotes pitch 
aS a subscript, denotes roll 


ratio of magnitude squared of platform [ to 
magnitude squared of single leg [Pf 


as a subscript, denotes single leg 
wave spectrum 


35 


as a subscript, denotes total platform 


non-dimensional quantity determined by leg spacing 
and wave frequency 


angle of wave incidence measured from roll axis 
angle of incidence of center of wave spectrum 


transfer function from incident wave amplitude to 
force (or moment) exerted upon the structure 


Dirac delta function 

wavelength of incident waves 

standard deviation of wave directional spectrum 
summation 


phase shift of the wave incident on leg i as compared 
to the wave at the center of the TLP 


normal distribution function 


radian frequency of incident waves 
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Appendix I. Computer Program Listing 


A. Discussion. 

Because of the complicated forms of the integrands of 
Equations III-17 and III-19 no closed form solution to the 
integral was found to exist. Yet the trigonometric func- 
tions with trigonometric arguments caused such rapid oscil- 
lations of the integrand that most numerical methods failed 
to converge. The method found to be the most useful was 
Gauss-Legendre quadrature (see Appendix III). A quadrature 
of fifth order was found to give sufficient numerical 
accuracy and is contained in the subroutine GAULEG. 

Because the wave directional spectrum was only non- 
zero between plus and minus 1/2 radians a correction factor 
was necessary to ensure its area was unity. This correction 
was accomplished by an eighth order polynomial approximation 
to the area under the normal curve. This approximation was 
developed by the author and is contained in the subroutine 
ECE: « 

Because of peculiarities in the computer system 
available the platform transfer functions are entered as 


function subprograms HEAV and FHEAVE for heave, and ROLL 
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and FROLL for roll, the former of each pair being the 
transfer function with uniform wave spreading. 
The program is written in FORTRAN-IV language and has 


numerous comment statements to aid in its understanding. 


B. Program listing. 


A listing of the programs for heave and roll, with 


associated subprograms, follows. 
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_— = a = _—— ts. ema TEE IE ee nn Se SE A A A REA a ANY -e Sp aetna cere cee renee ee 
Too ico ee 8 1 ORM 
esi LO 
en ea Cee gr AV 


Seti Ore) DSVIRTIO* (SD),LEG SPACING (D),AND SEPIOD CPER) 


ND=152. 
£9=2).55 
Pre=2. 


BereceLTc ViREANCE AND T 


VRI=5D"*2 

M=1.726533 "D/P rR**2 
Neos (5,5) 

eee eee te, SQUARE TEISION LES SLATEFORY ,AESVE EXCITATION”) 
wets €5,19) 

2O°48T C* *,515,°SRUSSTIAN WAVE DISTRIBOTINN’//) 
Sears (5,15) 0 

ee eto, LES SP8CT NGS", 749,F5.9,° FEET® ) 
mod. (S522) Pag 
ead (teers cf, * Pert 
Wes (6,75) 3D 
eee is, oFF SOAP D DEVESTIOV=°,1T29,F* .2,° RADIANS’//) 


‘—? 


epee sg ets) SGC ONO G ) 


Sor eeeulOne=Ncl 2 FO2 LOSS Sf TAILS SF  SRUSSIAN : 
H=7.579 796/85) 
Ehr ea nord ,P41t) 
PO2 ame = 1.) 
bat yo) S 
Gm oece ono c th ON ERETOP O05 TS LoSs OF TAILS=° ,F6.3) 
ero Ue SORE ADDING 


“DO 730909 25 
an) oi 


3 4-6 Ferd 


a =i t* Oo 


2 


-1.5797946 
0570796 


e _ 


2) 


SANTES CLL, XY,HTAV RIP ,27,VaP TT) 


Wit i wou 
i 
-) % 


saute 


Fema Gs 725) 213 

FORMAT (° *,TIS,*FACTO® DUE 7D WNIFORM SOREZDING (13)=",E11.4) 
WEIT (6,88) TP 

eer 2 >, (ERO CONTITIOV=",13,°)°/) 

dEITS (5,05) 


Meee et tg PETE T23,°51° ,rU5, “HN=t1/T3°,157,°SRROR?*/) 


229 
32) 
5 
1/8) 
15 
20 
25 
33 
55 
49 
& 


Mepeeer S erL OY OF CENTAR AT ANVE STSCTSOY 


Too. 


QAGhSts%, 

COzEIDES*1.735329E <2 

£O1=7 L490 

YVZS¥I+e+Ac 

Sera es (fb *, X01, Fae ay,Y,1EP,80,VAP,T) 
peav=y*c/S0 

UYSRTI/RI3 

Bem’ (6,5)) SOCES, 874, Fi, TEP 


Sete CC *° #8 TEI, tT S11648,TU2,FI10.43,7S2,13) 


ASNTISAINTS+S. 

Ee (h00Re .LR. 45.) 29 TO 199 
ST] 

END 
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~_A- ae 
we ree ~e —— ~ 


c TENSION LEG FLATFORWP 
C ROLL 
EXTERNAL ROLLeFROLL 
C SET STANDARD DEVIATION (SO) eLEG SPACING (DO) -¢4NGH FERTOD CPER?} 


D=1406 
SO=C e6S 
PERS2e 
C CALCULATE VARIANCE ANDO T 
2¢6c VAP SSDaed 
300 TOL e2269038*C/PER 222 
WEITE (€6¢5) 
= FORMAT C°1%eT1O4e*SGUCRE TENSION LEG PLATFORM, ROLL EXCITATION®) 
WRITE €560e19) 
ic FORMET €*% %,TISe*®GAUSSIAN WAVE OISTRIEUTION?S/) 


WRITE (6415) 2 
Us FOSMET C8 eT 1S5e*lLEG SFECING= "s1o0sFaeGy Seo e in 
WRITE (6925) PER 
2C FOR“MAT €® Se JlSe*PERTOC=%e TEC er Sele® SECCKHES*) 
WRITE (6025) SO 
ae FORMAT CC® eT1lSe*STENDARD CDEVIATIONHRN"%s 12S eFGe2e" RADIANS * 77a 


aQAn 


CORRCCTION FACTOR FOR LOSS OF TAILS OF SAYSSTIAN 
WINTH=S1-570796/5SN 
CALL ERF C4IDTHe PHI) 
CrLeFC2erPHJ m1.) 
WRITE (6¢e36) C 
30 FORMAT C€*® *eTISe*COPRECTICON FACTOR CUE TO LOSS GF TAILS=%eF G23) 
FACTOR SUE TO UNTFORM SPRE ACING 


ry 


XL="-1-e5739796 
XU= 1.570796 
ecsde 
call qaulea Cxlexusrollereire=CeVAPReT) 
RIZ=R*eDex? ’ 
WRITE (6e35) RI3 

25 FORMAT €% %.725¢ °F ACTOR CUE TO UNIFORM SPREADING CI3)=% £11.49) 
WRITE €6e4f) Te ’ 

46 FORMAT €® %eT20e*(CERPOR CONDITION=*-eI3/) 
WRITE (6445) 

45 FORMAT €® %9% BLTA% se TAZ oF T1*%* eg F459 *HNSETIST 2% TE2 0%“ RROR? */) 

C 
C VERY GIFPECTION OF CENTER OF WAVE SPECTRUM 

BCOEG=O0e6¢ 

100 PO=8OCEC#12745329E-2 
XLI=XL*39 
XU1=XUeB8O 
call caulead (xtlexulefroll eyeiere ite VAR oT) 
RIL1IsY¥*«PewoeC/EO 
HRERTIZRI3 
WRITF (6050) BNDEGeRITIeHNeIER 

SO FOPMAT € %eF5eleTI7eELeSeTA2eF LD Ge TESe 13) 
FANEGSPODEGSE. 
IF (€8CDEG eLEs S00) GG TO J00 
STGP 
END 


a 


@ |) 
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SUS OO Ge ILS CVE te Cre Ye LEP FN ,VEPe.T) 
See Pe ORDER SGAUSGS@-LESEIDPE QUADTURTURE 

T=0 

H=(XN-XL)/15.- 

A\=eee 2 sci le 

R2=9.47862797 

R3=0, 2628849 


X1=5.9041798 
Mm2Q=7.529UF93 


SU“ =C. 


wo ,57 

Y2=.5*H 

A=XL4"%2 

Y1I=X1*F2 

Y2= 2°02 

T=I+1 

So" 1=suU"™ 

SU¥s% © 

P1="+¥4 

RB2=4-Y1 

R3=A+42 

PY=h-V2 

SUM=SSUMFAT* CTITC81,30,7AP,TIFPCTOR? 66 ,VAR,T))+h2* 
CHeGORS Zoo gyi | (SP eTCRU RO TAC GT) yeas*FCTC"® ,RO,7AP TT) 


A=A +4 

IFC XUMA 4,0 ,3 

SUr=e oS IM 

PME ) . 36) GQ TOF 
ERP=ARS((SU%=SI%1) /SU™) 


IF (1.67.10) 3 
Bee CPR R= sc) Ss 


Ol ae alte) 
rao a 
TERS4 

y=si™ 

RETUQY 

TER=1 

Y=SU™ 

RETURY 

ED 
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| 
+2 
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Ds 

lay oe ee 
PRi=%. 

Cae ae 

ee 

PUT= 5 

NSA (Giese) 20g) Sy 
MTN CON = 1 a a 
NTY (3) =-6.1512795=2 


IV (b )e= 1.732 es —2 
\UY Ce j= > See 
ie 


3 
RIV (T)= 1,4392515-3 
5 


BINViCR )=-9 6915 29r— 

AX tf I=7,2 

Va *y 

DHT =rAT eve aT? (7) 

Ta oe eee 
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Meee LON FHeadV (368 Ge VAP od 

eM te ONPG He (CP eo tC TOP Tt CAG CRS Pde Ch et EME tT eS OT) dD)» 
eee (NK ee OS (C2 et VE")) 

RETURN 

ED 


Mere O3On HEAVY (bo PG eVAts 7) 
MEE 75246 1.4 COS 6 TeCOS (RF) )* 0104015 (74ST 138? 
meen E NN 
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Seer i ON T8cuL(r,286,V3P,7) 
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7 S720 -( 2-27) **2/(9,*V2P)) 
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ap 


Spee TTC) 297L(2,2°,¥2°,T) 
RIL L= 0 3183 ID*( THOM (THOT (8)) HOT -TOSETHS E(B) )) 


2eTye" 


ED 
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Appendix II. Directions For Use of Program 


A. “Discussion. 

This program was intentionally designed to remain as 
much as possible in discrete sections, each performing a 
specific purpose. In this manner the program can be more 
easily understood by others and converted to use for many 
different sizes and types of structures. 

An error flag has been added to the subprogram GAULEG. 
If GAULEG is unable in a reasonable time to get consecutive 
iterations to agree within one percent it flags this fact 
to the user by indicating an error condition of one. If the 
error condition is zero then GAULAG has gotten the same an- 
swer (within one percent) for two iterations with the 


second iteration having twice as many increments as the first. 


B. Numerical parameters. 


To make it easier to understand the program flow the 
leg spacing (D), standard deviation (SD), and period (PER) 
are simply assigned values in FORTRAN statements. They could 
just as easily be read from data cards, input from an inter- 


active terminal, or varied in do-loops, whichever is most 
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convenient for the user on the available system. 

To adopt the program to a square TLP of interest one 
must simply replace the assignment statements for leg spacing 
and for the period of interest. If many periods are in pOmee 
ant then a do-loop for period should be used returning to 
statement 300 to ensure the value of T is updated. If a 
do-loop is used to vary standard deviation it must return 
to statement 200 to revise the value of the varience. 

The center of the wave spectrum (3,) was varied in an 
implied do-loop with statement 100 because one is normally 
interested in the response to waves coming from all directions. 
Because of the symmetry properties of a square TLP the angle 
must only be varied over 45 degrees for heave and 90 degrees 


BOmerOll. 


C. Conversion to other types of platforms. 


The information concerning the platform geometry is 
contained in the platform transfer functions. In the computer 
program the definition of the functions in the function 
subprograms (ROLL, FROLL, HEAV, and FHEAV) entirely defines 
the structure. Note that a factor of d* was left in the main 
program rather than ROLL and FROLL since this factor will be 


present in moment transfer functions for any type of structure. 
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Thus to change this program to accomodate any shape TLP 
one must simply perform derivations similar to those in sec- 
tion III-A and use the results to replace the function defin- 


itions in the function subprograms. 


D. Conversion to other wave spectra. 


This program was explicitly written for a normally dis- 
tributed wave directional spectrum because it is so easily 
varied. Simply by varying the standard deviation one can 
obtain as directional or uniform a spreading as desired. 
However, if another spectrum is preferred this program can 
be adapted. 

The most obvious place to change is in the definition 
of FROLL and FHEAV. There the new spectrum would be substi- 
tuted for the exponential term. In the main program the 
correction factor C/SD is applied after calling GAULEG using 
FROLL or FHEAV. This would be changed or eliminated. Since 
C, the factor for loss of tails, is no longer necessary the 
entire subroutine ERF could then be eliminated. New correc- 
tion factors for the new spectrum would then be introduced. 

It is also possible to normalize the transfer function 
using something other than uniform spreading. For this ROLL 


and HEAV only must be revised. If correction factors are 
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necessary they are best introduced in the main program after 


calling GAULEG for the integration. 
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Appendix III. Gauss-Legendre Quadrature 


A. Legendre Polynomials 


1. Definition. Define a set of polynomials P,(x) of 


degree n, where n is a non-negative integer, such that 


: 
ft Pim (X) Pf) dx = cy) on 


and Poh) el, 


The first is an orthogonality condition and the second 
is necessary to fix c,. Together they define a unique set of 
polynomials known as the Legendre polynomials, the first few 


of which are 


Py 6x) =z ] 

P, (x) se 

Po(x) = 4(3x?-1) 
P, (x) = &(5x'-3x). 


In general they can be found by 


2 qn 2 n 
Se) cee (x*—e1l)™, 


P(x) 


Note that Pfs) has n zeroes in the interval from minus one 


to one and is even for n even and odd for n odd. 
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2. Polynomial approximation. Any function piecewise 


continuous in the interval from minus one to one can be 
expanded in terms of an infinite series of Legendre polynomials. 
It can be approximated by a finite series of Legendre poly- 


nomials. 


co n 
y (x) = % Sm P,, (x) = eg CM Pi, (x) 


B. Integral approximation 


A simple means to evaluate the integral in the region 


minus one to one is as follows: 


i 
J y(x) dx = f Cm Pm(x) dx 


omg 


Multiply both sides by 1 = P)(x): 


i y(x) dx = r Cc 
1 0 


tnt Cy one 


PB. P,, (x) P, (x) ax 


= Co f P(x) P(x) dx 
i y(x) dx = 2c 


Thus the problem has now become one of finding Cau the 
coefficient of P(x). One possible method will be shown 


using the example where n=2. 


We desire to find coefficients and arguments such that 
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| V(x) dx@= Cae A, y (x) ) * A, y(x,) 


For any two points. x, and x,, we expect 


Cee exe GC EF (x) ) eevee Sl = yA) 
Cc, OP, (x, ) aa Cc, P. (x) + C, Ee (x) = y(x,) 
Noting that P(x) = 1 and P, (x) = x and choosing x, and 


x, to be the zeroes of P, (x) this becomes 


Q 
oe 
Q 
ra 
il 


y (x) 


Q 
+ 
Q 
a” 
i 


y (x) 


Where x, and x, are + (1/73) . By adding these two 


equations and noting that x, = -x, we have 


2C. = y(-1//3) + y(1//3) 


0 


which is our desired answer. 
To get a larger degree of accuracy a larger n 1s required. 
This method for determining 2 C, remains effective. For the 


case of n = 5 the solution is (to seven significant digits) ? 


1 
Fy (x) ax 


(0.568 888 9) {y(0)} 


op 


(0.478 628 7) {y(.538 469 3)+y(-.538 469 3)} 


+ (0.236 926 9) {y(.906 179 8)+y(-.906 179 8)} . 
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However, one is not always interested in the interval 
from minus one to one. A transform must then be used to 
transfer the region of interest to that interval. For addi- 
tional accuracy without the need for a larger order quadra- 
ture the total interval can be divided into subintervals. 
Each subinterval can be transformed seperately and the 
result summed. Thus each subinterval will be approximated 
by a fifth order polynomial. This is the method employed 


in the subroutine GAULEG. 
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